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Abstract 

We have studied the existence of self-dual solitonic solutions in a generalization of the Abelian 
Chern-Simons-Higgs model. Such a generalization introduces two different nonnegative functions, 
wi(|0|) and a;(|0|), which split the kinetic term of the Higgs field - |D M </>| 2 —> oii(|(/)|)|Z)o</>| 2 — 
Lo{\4>\)\Dk(j>^ - breaking explicitly the Lorentz covariance. We have shown that a clean implementa¬ 
tion of the Bogomolnyi procedure only can be implemented whether w(|0|) oc /3|</>| 2/3_2 with /3 > 1. 

The self-dual or Bogomolnyi equations produce an infinity number of soliton solutions by choosing 
conveniently the generalizing function wi(|0|) which must be able to provide a finite magnetic field. 

Also, we have shown that by properly choosing the generalizing functions it is possible to reproduce 
the Bogomolnyi equations of the Abelian Maxwell-Higgs and Chern-Simons-Higgs models. Finally, 
some new self-dual \<j >\ b -vortex solutions have been analyzed both from theoretical and numerical 
point of view. 

PACS numbers: ll.10.Kk, ll.10.Lm 


I. INTRODUCTION 

A time ago it was shown that (l+2)-dimensional mat¬ 
ter field interacting with gauge fields whose dynamics 
is governed by a Chern-Simons term support soliton so¬ 
lutions [1], [2] (for a review see Refs. [3], [4], [5], [6] 
and [7]). These models have the particularity to become 
self-dual when the self-interactions are suitably chosen 
[8], [9], [10], [11]. When self-duality occurs the model 
presents interesting mathematical and physical proper¬ 
ties, such as the second order Euler-Lagrange equations 
can be solved by a set of first-order differential equa¬ 
tions [12], [13] and, the model admits a supersymmetric 
extension [14]. The Chern-Simons gauge field dynamic 
remains the same when coupled with matter-fields either 
relativistic [8], [9] or nonrelativistic [10], [11]. In addi¬ 
tion the nature of the soliton solutions can be topological 
and/or nontopological [15]. 

The inclusion of non-linear terms to the kinetic part 
of the Lagrangian has interesting consequences, as for 
example, the existence of topological defects without a 
symmetry-breaking potential term [16]. In the recent 
years, theories with nonstandard kinetic term, named k- 
ficld models, have received much attention. The /e-field 
models are mainly in connection with effective cosmolog¬ 
ical models [17], [18], [19], [20], [21], [22], [23], as well as 
the tachyon matter [24] and the ghost condensates [25], 
[26], [27], [28], [29]. The strong gravitational waves [30] 
and dark matter [31], are also examples of non-canonical 
fields in cosmology. The investigations concerning to the 
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topological structure of the fc-field theories have shown 
that they support topological soliton solutions both in 
pure matter models as in gauged field models [32], [33], 
[34],[35], [36], [37], [38], [39], [40],[41], [42], [43], [44], [45], 
[46],[47], [48], [49]. These solitons have certain features 
which are not necessarily shared with those of the stan¬ 
dard models [50], [51], [52], 

The aim of this manuscript is to study a Chern-Simons- 
Higgs model with a generalized dynamics which breaks 
Lorentz covariance, i.e, 

\d„<p \ 2 wl (i<(>i)iA)/>r - u(\<t>\)\D i( t>\ 2 . a) 

The nonstandard dynamics is introduced by the func¬ 
tions ui\ and u>, depending on the Higgs field. During 
the implementation of Bogomolnyi trick is demonstrated 
that self-dual configurations exist if the function u> is pro¬ 
portional to \(j)\ 2 P- 2 with (3 > 1. On the other hand, the 
function uq remains arbitrary but near the origin should 
behave as \(f>\ 2S with 6 > —1 in order to have a well behav¬ 
ior for the magnetic field. In particular we have chosen 
the functions ui\ and ui, to be 

^(101) = |^| 2M , W (I0|) = (N + 1)|0| 2JV , (2) 

where M > — 1 and N > 0. This way, the Bogomolnyi 
equations produce an infinite number of soliton solutions, 
one for each value of the pair ( N,M ). It is possible to 
show that for particular values of A, M, the Bogomolnyi 
equations of the Maxwell-Higgs or Chern-Simons Higgs 
models can be recuperated. Finally, we have constructed, 
analytically and numerically, novel soliton solutions for 
some values of N and M. 
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II. THE THEORETICAL FRAMEWORK 


Following the same ideas introduced in Refs. [32], [33], 
[34], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44], 
[45], [46], [47], [48], [49], we start by considering a gen¬ 
eralized (2 + l)-dimensional Chern-Simons-Higgs (CSH) 
model where the complex scalar field possess a modified 
dynamic. Such a model is described by the following 
action, 


S = S cs + d 3 x 


- cornel 2 -vm 


( 3 ) 

where S cs represents the Chern-Simons action given by 


S cs = d 3 x-e^ p A p F„ p . 


( 4 ) 


however such as it happens in usual CSH model, the elec¬ 
tric charge is nonnull and proportional to the magnetic 
flux: 

Q = k Jd 2 x B = K<f>. (10) 

Therefore, independently the functional form of the gen¬ 
eralizing functions wi(|<^|) and u}(\4>\), the solutions al¬ 
ways will be electrically charged. 

Likewise, the Ampere law reads 

^F a0 + euJ k = 0. (11) 

Along the remain of the manuscript, we are inter¬ 
ested in time-independent soliton solutions that ensure 
the finiteness of the action (3). These are the stationary 
points of the energy which for the static 


The covariant derivative D p (f> is defined by 

D p <j> = d p (j) - ieA p (f>, (5) 

with p, = 0,1, 2. The metric tensor is g pv = (1, —1, —1) 
and e^ p (e 012 = 1) is the totally antisymmetric Levi- 
Civita tensor. 

In action (3) we notice the usual Higgs kinetic term, 
| D p (j)\ 2 = \D 0 (f>\ 2 — \D k (j>\ 2 , was replaced by a more gener¬ 
alized term, wi(|<(>|)|-Do<y!>| 2 — w(|((>|)|.Dfc<(>| 2 , which breaks 
explicitly the Lorentz covariance. The dimensionless 
functions uq(|<(>|) and uj(\(f>\) are nonnegative and, in prin¬ 
ciple, arbitrary functions of the complex scalar field <fi. 
The function F(|<(>|) is a self-interacting scalar potential. 

The gauge held equation obtained from the action (3) 
is given by 


1^F ap - ej p = 0, (6) 

with = ujiSq J° + ut6%J k is the conserved current of 
the model and J p = ^(D 1 *(/))*— cf>* (D 1 *<j))\ is the conven¬ 
tional current density. Similarly, the equation of motion 
of the Higgs held is 

0 = (doOJ^Dot + ojMDocj))-^Docj)] 2 (7) 

Blj dV 

~{d k u)D k <t> - ojD k (D k (j)) + — \D k c />\ 2 + 

From Eq. (6), the Gauss law reads 

kB = ecuiJo, ( 8 ) 

we observe the Gauss law of Chern-Simons dynamics is 
modihed by the function wi(|0|) such that now the con¬ 
served charge associated with the U{ 1) global symmetry 
is given by 


E = fd 2 x[—nA 0 B - e 2 wiA 2 |</>| 2 + u\D^\ 2 + V(\<f>\)] . 

( 12 ) 


From the statitic Gauss law, we obtain the relation 

k B 

An = - 


2e wi|<^| 


2 ’ 


(13) 


which substituted in Eq. (12) leads to the following ex¬ 
pression for the energy: 


E= d 2 


tv B 2 


■ uj\Di(j)\ 2 + V{\cj>\) 


(14) 


j |_4 e 2 wi |<(>| 2 
To proceed, we need the fundamental identity 

\Di4>\ 2 = \D ± ^\ 2 ± eB\(f>\ 2 ± ^e lk d t J k (15) 

where D±<p = D i</>± iD 2 </>. Then, by using (15), we may 
rewrite the energy (14) as 


E = d 2 


k 2 B 2 
4e 2 uji\cp\ 2 


vm+uiD^i 2 


2 1 

±eu>B \cj>\ ±-Loe ik d z J k 


(16) 


We observe that the function cc(|^|) in the term 
wei k diJ k preclude us to implement the BPS procedure, 
i.e, the integrand must be expressed like a sum of squared 
terms plus a total derivative plus a term proportional to 
the magnetic field. Therefore, the key question is about 
the functional form of w(|^!>|) allowing a well defined im¬ 
plementation of the BPS formalism. We start the search¬ 
ing of the function w(|</>|) from the following expression: 

^ik'Bi k ) — M€i k BiJ k T £i k {BiUJ^J k . (L7) 
By manipulating the last term ei k (diU})J k it reads 


Q 



x eu)\J 


o 


( 9 ) 


£ik k 


But 


m 


2 (pi l^l 2 ) dkJk, 


(18) 
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where we have used the fact of w be a explicit function of 
\(j)\ 2 . After some algebra the term e ik (d/w) J k becomes 

Buj Blj 

e ik {diLj)Jk = |0 | 2 —^ e ik d l J k + 2 eB |0 | 4 ——(19) 

d 101 <9 101 

Substituting this equation in (17) we arrive to 

( duj \ Ouj 

+ l ^| 2 Q\(j>\ 2 ) 2 el ?|(/)| 4 ^| 0| 2 • 

( 20 ) 

Here we impose that the function ui satisfies the fol¬ 
lowing equation: 


W+ |0| 2 7^2 = 




( 21 ) 


with p a real constant. By solving Eq. (21) we obtain 
the explicit functional form of u/(| 0 |), 


w = c]^- 2 , 


( 22 ) 


where the constant C adjusts conveniently the mass di¬ 
mension of OJ. 

The key condition (21) allows to rewrite Eq. (20) in a 
more simplified form 

e ik di(ujJk) = Puei k diJ k + 2e(/3 - l)uiB\4>\ 2 , (23) 

allowing to write the term uietkdiJk in the following way 

ue ik diJ k = ^ e ik di (wJ k ) - 2 \<t >\ 2 . (24) 

By introducing it in Eq. (16), the energy becomes 
k 2 B 2 


E = d x 


4e 2 u/i|0 | 2 


Vm+u\D±<t>\ 2 (25) 


±^euB |0| 2 ± e ik di (wJ fc ) 

P ^P 


We write the two first terms as 

k 2 B 2 k 2 1 

+ P = 


4e 2 |0| 2 u/i 


4e 2 |0| 2 wi 

/T 

e| 0 | V wi 

By substituting in (25), we have 


Bt ^ |0| 


(26) 


E = d 2 x 


4e 2 |0| 2 o / 1 ( ST 

+w|£> ± 0| 2 ± ,/. <9, (wJfc) 

F 1 


±s * -rn\/—+ -p euJ W 
e| 0 | V wi p 


equal to ev 2 , it allows to define explicitly the form of the 
potential H(| 0 |), 


4 4 

e u 


^( 101 ) = - 5-^101 


1 - 


101 


2/3' 


i; 2 /3 


(28) 


where we have substituted the explicit form of w(] 0 |) 
given by Eq. (22) with C = pv 2 ~ 2 ^ in order to the 
vacuum expectation value of the Higgs field to be |0| = v. 
The function uq (|0|) still remains arbitrary. Hence, the 
energy (27) reads 

E =jd 2 X {±e V 2 B±^ k d^Jk) + .\D^\ 2 (29) 

2" 


4e 2 |0| 2 cji 


9p^7; 2 

H T ^^i|0| 2 


1 - 


101 


2/3' 


i; 2 / 3 


We see that under appropriated boundary conditions the 
total derivative gives null contribution to the energy. 
Then, the energy is bounded below by a multiple of the 
magnetic flux magnitude (for positive flux we choose the 
upper signs, and for negative flux we choose the lower 
signs): 


E > ±ev 2 / d 2 xB = ev 2 \<&\. 


(30) 


This bound is saturated by fields satisfying the Bogomol- 
nyi or self-dual equations [ 12 ] 


£>±0 = 0, 


B = ± 


2e 3 v 2 


-oj 1 


1 - 


101 


2/3' 


,;2/3 


(31) 


(32) 


If we require that the magnetic field be nonsingular at 
origin, the function o/i(| 0 |) should behave like | 0| 25 with 
8 > — 1. On the other hand, positivity and finiteness of 
the BPS energy density requires /3 > 1. 

Below we study interesting models by given a specific 
form of the functions uj and ui 


III. SOME SIMPLE MODELS 

In the following we analyze some interesting but simple 
models by setting, 


a/(|0|) = (N + 1). 


,|2 N 


„2 N 


U/i(|0|) = 


101 


2 M 


„2 M 


(27) The BPS potential (28) reads 


To finish the BPS procedure, we observe that if in the y (| 0 |) _ 

third row the term multiplying to the magnetic field is 


4„6 UI2M+2 


e 4 u 6 | 0 | 


c 2 ni 2 A ^-)-2 


i - 


101 


2/V+2 ' 


I.2/V+2 


(33) 


(34) 
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and the BPS equation (32) becomes 


B = ± 


2e 3 v 4 \(j)\ 2M+2 
1^2 -\-2 


\<fi\ 2N+2 \ 

v 2N+2 J ■ 


(35) 


Here, it interesting to note that for TV = 0 and M = 0 
the self-duality equations (31) and (35), becomes the well 
known Bogomolnyi equations of the Chern-Simons-Higgs 
theory [8], [9], 

9 p 3 

D±4>= 0, B = ± — \<P\ 2 {v 2 - \</>\ 2 ) ■ (36) 

Av 

In the case N = 0 and M = — 1, we have 

D±4> = 0, B = ±=^(v 2 - \0\ 2 ). (37) 

These equations are essentially the Bogomolnyi equa¬ 
tions of the Maxwell-Higgs model, whose solutions are 
the well known Nielsen-Olesen vortices [53]. The differ¬ 
ence lies in the fact that, here, our self-dual solitons not 
only carry magnetic flux, as in the Higgs model, but also 
U{ 1) charge. This is a consequence that in our theory the 
dynamics of gauge field is dictated by a Chern- Simons 
term instead of a Maxwell term as in Maxwell-Higgs 
theory. So, for N = 0 and M = —1, we obtain self¬ 
dual configurations which are mathematically identical 
to the Nielsen-Olesen ones but differently our solutions 
have electric charge. 


A. Vortex configurations 


Specifically, we look for axially symmetric solutions us¬ 
ing the standard static vortex Ansatz 

<t> = vg(r)e ine , = - ° (r) ~ n (38) 

er 

The Ansatz allows to express the magnetic field as 

B = ~— (39) 

er 


where ' denotes a derivative in relation to the coordinate 
r. Likewise, the BPS equations (31) and (35) are written 
as 



(40) 


B =--= ±'^-^g 2M + 2 (l - g 2N + 2 ) . (41) 

ern z 

These equations are solved considering the profiles g 
and a are well behaved functions satisfying the follow¬ 
ing boundary conditions 

3(0) = 0, «(0) = n, (42) 

(43) 


The BPS energy density of the model reading from 


GBPS 2" 


dr re F 


(44) 


is given by 


_ 2e 4 V 6 2M+2 (-> 2N+2\ 2 

^bps 2 y v y j 

hj 

+2v 2 (N + 1) g 2N , (45) 

the requirement of finite energy density, for all values of 
the winding number n, imposes N > 0 and M > —1. 


B. Checking the boundary conditions 


We obtain the behavior of the solutions of Eqs. (40) 
and (41) in the neighborhood of r —> 0 using power series 
method, 


g(r) = G n r n 


P 4 „ 4 (n \ 2M+3 r n(2M+3)+2 

6 W {Un) ---.-+ ■ • ■ (46) 

2 k 2 (nM + n + 1) 


a(r) 


e 4 v 4 (G n ) 2M+2 r ra ( 2M + 2 )+ 2 
k 2 (nM + n + 1) 


(47) 


It verifies the boundary conditions given in Eq. (42). 

For r -A + 00 , the behavior of the soliton solutions 
becomes similar to the Nielsen-Olesen vortices, 


P Ill'S' 

~ ‘- C - -Jr 

(48) 

a(r) ~ G ao m s ^Fe- m * r , 

(49) 


where G ^ is a numerical constant determined numeri¬ 
cally and m s , the self-dual mass of the bosonic fields, is 
given by 

m s = V-ZV + l. (50) 


It is verified that for N = 0, the mass scale is exactly the 
one of the Chern-Simons-Higgs model. 


vm 


e 4 v 6 \(/)\ 2M+2 ( i ^ r + n 2 

k 2 y 2M+2 I 1 y 2N+2 J ’ 


(51) 


C. Numerical analysis 

Below, without loss of generality we set e = 1, v = 1, 

K = 1. 

Before performing the numerical solution of the self¬ 
dual equations (40) and (41) we do the following obser¬ 
vations in relation to the BPS potential (34): First, it 
provides a \(j)\ 4 potential for M = —1 and N = 0, 

nw) = (i-w 2 ) 2 , 


g(oo) = 1, a(oo) = 0. 


(52) 
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Second, the BPS potential also provides a family of 
\(f>\ 6 potentials when the condition M = —2TV is satisfied 
and N is restricted to the interval 0 < N < 1/2, 

m^i) = H" 4 JV+ 2 (i-H 2JV+2 ) 2 - (53) 

Below, our numerical analysis considers only these two 
potentials to solve the BPS equations (40) and (41). In 
particular, we solve the Bogomolnyi equations only for 
winding number n = 1. 

In the figures, the red line represents the case M = — 1 
and N = 0 providing the Nielsen-Olesen-like vortices, 
whereas the blue lines depict the vortex solutions for the 
values of M and N generating some |0| 6 -potentials. In 
particular, we have plotted three solutions in blue lines: 

• M = 0 and N = 0, which generates the well know 
Chern-Simons-Higgs vortices. 

• M = —0.5 and N = 0.25, associated to the self¬ 
dual potential 

v(\4\) = \<f>\ (i- I#) 2 - (54) 

• M = — 1 and N = 0.5, associated to the self-dual 
potential 

vm= (l-|^| 3 ) 2 . (55) 

Note that in the cases M = — 1, N = 0 and M = — 1, 
N = 0.5, i.e. the cases associated to the potentials (52) 
and (55), there is only one degenerate vacua at \<j>\ = 
1. This fact leads us to similar solutions, which can be 
appreciated in Figs. 1, 2, 4, 3. 

For the cases where M ^ — 1, the \(f>\ 6 potential have 
two vacua: \<j>\ = 0 and \<j>\ = 1. In these cases, the 
profiles of the magnetic field are rings whose maximum 
amplitude, for increasing values of N, approaches to the 
origin (see Fig. 4). Also the profiles of the BPS energy 
density have a ring-like format ((see Fig. 5) and the ring 
format is explicit for n > 1. 

On the other hand, for electric field, whenever the val¬ 
ues of M and N here considered, the profiles always are 
rings around the origin (see Fig. 3). 

IV. REMARKS AND CONCLUSIONS 

In summary, we have proposed a generalized abelian 
Chern-Simons-Higgs model with explicit breaking of 
Lorentz covariance and explored the respective Bogomol¬ 
nyi framework. During the implementation of the BPS 
trick it is shown that the generalized functions should 
satisfy some requirements: The function w(|<^>|) must be 
a monomial, i.e., w = C\<p\ 2/3 ~ 2 for all P > 1 and the func¬ 
tion uji(\(f>\) must be regular at the origin (uq oc \(f>\ 2S with 
5 > — 1). Under such conditions imposed on the general¬ 
ized functions, it is guaranteed the existence of self-dual 



r 


FIG. 1: The profiles of the Higgs field g(r) for n = 1. The 
red lines represent the solutions for a |0| 4 potential and blue 
lines for \(j >\ 6 potentials. 



r 

FIG. 2: The profiles of the gauge field a(r) for n = 1. The 
red lines represent the solutions for a \(j>\ A potential and blue 
lines for |(/>| 6 potentials. 



r 


FIG. 3: The profiles of the electric field E(r) = A' 0 (r) for 
n — 1. The red lines represent the solutions for a \(j>\ A potential 
and blue lines for |0| 6 potentials. 
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FIG. 4: The profiles of the magnetic field B(r) for n = 1. 
The red lines represent the solutions for a |</>| 4 potential and 
blue lines for \4>\ 6 potentials. 



FIG. 5: The profiles of the BPS energy density e bps (r) for 
n — 1. The red lines represent the solutions for a |</>| 4 potential 
and blue lines for \(p\ 6 potentials. 


solitonic configurations satisfying Bogomolnyi equations 
whose magnetic field and BPS energy density are well 
behaved. As we expected, the infinity family of self-dual 
configurations have finite energy which is proportional 
to the magnitude of the magnetic flux. In particular, we 
have studied the self-dual vortices provided by the choice 
wi(|</>|) = (N + 1) \4>\ 2N and u>i(|0|) = |0| 2M . It was 
shown the vortex solutions of the Maxwell-Higgs model 
and the Chern-Simons-Higgs model can be also obtained. 
Besides that, we have constructed two new solitonic solu¬ 
tions which correspond to Chern-Simons theory coupled 
to two types of |0| 1 2 3 4 5 6 7 8 9 10 11 potentials given by Eqs. (54) and 
(55), respectively. 

Finally, it is worthwhile to point out that existence 
of BPS states is linked to the existence of a M = 2- 
extended supersymmetric model [54], We are studying 
such a possibility despite the fact that in this model the 
Lorentz symmetry is explicitly broken. Advances in this 
direction will be reported elsewhere. 
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